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1 Introduction 

In this work we analyse for the first time the effect of the infrared structure of QED on the entanglement of charged 
qubits. 

The qubits used in Quantum Information Theory are built from selected degrees of freedom which belong to (more 
or less complex) quantum systems. One of the simplest cases could be the spin of a charged particle or the photons 
polarisation. A new discipline is born in order to properly analyse these systems, it is what Peres and Terno [1J called 
Relativistic Quantum Information Theory, which concerns the role of relativistic-like considerations in the Quantum 
Information Theory. 

The attempts to include the dynamical effect of Quantum Electrodynamics (QED) on the correlations between 
particles are very recent. We might consider that they started with the studies of the effect of the QED spin-spin 
interactions on the entanglement [2 and the violation of Bell Inequalities due to QED 3 . Complementarily, an "infi- 
nite" entanglement generation between interacting systems momenta has been predicted [4], as well as entanglement 
transfer between the space degrees of freedom and the qubits. Both cases in the QED frame. 

However, the formulation of Quantum Information Theory in the relativistic context is still far from being com- 
plete. For instance, there is no satisfactory notion of localization [5], [B], [7J or, for example, it is not clear what would 
be the proper spin operator in order to characterize that degree of freedom in the context of Relativistic Quantum 
Information Theory [8]. 

In the last decade, starting from the work of Czachor [3] and, afterwards, a lot of papers such as [I], [TU], [JJ, [TT] . 
[12j . |13j . |14j.|15j and many others have considered relativistic effects on Quantum Information Theory, proposing a 
new phenomenology regarding the entanglement from a Quantum Field Theory approach. 

One of the main problems of this kind of works is the need to use perturbation theory so as to characterize 
interacting and space-like separated qubits. The simplest way to do it is to describe them by means of free asymptotic 
states ( "in" and "out" ) as it is common in Scattering Theory, postulating that in the large time limit the interaction 
is adiabatically deactivated. With this assumption, the free field operators in the interaction picture would give us 
the description of physical particle states that properly evolve in the asymptotic regime. 

Nevertheless, it is well known [IBJ that the free Hamiltonian is not the large t limit of the QED Hamiltonian, 
therefore it does not make sense to identify the free states b\(p) |0) (where b\(p) is in the interaction picture) 
with asymptotic states within QED. Complementarily, divergences appear in every order of the perturbation theory 
because of emission and absorption of soft photons (photons whose energy and momentum are much lesser than the 
characteristic masses and energies of the processes) . However, the divergences that come from real soft photons are 
canceled by the divergences that come from virtual soft photons so that the cross sections are finite [TT] . 

On the other hand, the states that evolve properly in the asymptotic regime are safe from infrared divergences 
16J although, as it will be shown, the evolution under the electromagnetic interaction in the asymptotic limit makes 
it very difficult to identify charged particle states, actually, an undetermined number of soft photons appear, spoiling 
the gauge-invariance of the asymptotic states. What is more, in some papers it is said that the very concept of 
charged physical particle may not make sense rigourously speaking. This fact would be disastrous when you are 
trying to build qubits from charged particles. 

This fact is directly related with the masslessness of the photon, it leads to a series of pathologies in QED whose 
effect on the qubits and their entanglement is still unknown: isolated poles (associated with the mass of the charged 
particles) do not appear in the energy spectrum of QED; as a matter of fact, corresponding to each one, there appears 
a threshold of continuous states corresponding to states with an arbitrary number of soft photons. In those conditions, 
the electron propagator no longer corresponds to a simple particle pole (p 2 — to 2 ) -1 , instead, the exponent acquires 
an irrational value (which depends on the fine structure constant) that leads to a cut in the complex plane (concept 
of infraparticles |18j.) 

Along the second half of the XX century it has been discussed whether physical states could be associated to 
charged particles or not, and, if they could, what would be their nature. The first contributions to this topic come 
from Dirac himself [19 . 

The effect of the infrared structure of QED on the entanglement of a system of two charged particles (and which , 
therefore, undergo electromagnetic interaction) has not been considered yet, although Peres and Terno [TJ identified 
it as a problem, noticing that the concept of charged qubit is merely an idealization, and then, the properties of an 
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entangled state of charged particles would be modified in an unknown way due to the presence of an undeterminate 
number of soft photons. 

Recent works [20] , [21] have gone further into the procedure called "dressing" : construction of asymptotic states 
with good enough properties to be identified with physical charged particles by means of restoring their gauge 
invariance. However, no one of the previous works have treated in depth the phase effects and divergences which 
appear derived from the asymptotic evolution |16j . [T] when the asymptotic states are built. In this work these topics 
are dealt with. 

The objective of this work is, therefore, to characterize the effect of soft photons on the entanglement of two 
charged qubits, something that has not been done yet. In order to achieve it, we will need to extend the theory 
associated with the construction of asymptotic states and the "dressing" of these states so as to apply those results 
to our case. Here, we will compute for the first time the effect of the infrared structure of QED on the entanglement. 
Eventually, we will prove that all the correlation functions of the spin degrees of freedom (and particularly, EPR 
correlations) are not modified at all by the presence of soft photons when we are considering charged qubits built 
from proper asymptotic particle states. 

In section 2, we raise the most general two charged particles state in the free case. In section 3 all the aspects 
concerning asymptotic dynamics (and the construction of asymptotic states) are treated. In section 4 we apply the 
"dressing" procedure in order to build physical states (gauge invariant) from asymptotic charged particles states. 
Eventually, in section 5 the entanglement of these physical asymptotic states is compared with the entanglement in 
the free case. Conclusions are summarized in section 6. 



2 Two Dirac's Qubits. Free case 

Let us start with the description of two charged qubits. We are only concerned with their spin degrees of freedom. 
As charged particle we will understand (in this section) the idealized concept of free particle [T] which is commonly 
used in the frame of Rclativistic Quantum Information Theory. 

The most general state of 2 free fermionic qubits may be written a^] 

I*) = / d3 Pl d3 P2 ( P<J 1 a 2 (Pl,P2) |(Ti,pi ; CT 2 ,P2> (2.1) 

Where we have defined 

\u uVx ■ a 2 ,p 2 ) ee \ax,a 2 ) \ Pl ,p 2 ) ee &£> 2 )&i>i) |0) (2.2) 
Being 6j.(p) a fermionic operatoij^] in the interaction picture. It must be satisfied that 



J d 3 p 1 d 3 p 2 \ip ailJ2 (pi,p 2 )\ 2 < oo 



Let us build the density matrix associated with this general state 

|*)(*| = J d 3 p 1 d 3 p 2 d 3 p[d 3 p 2 S ai „ 2t „ [ „,(p 1 ,p 2 ,p' 1 ,p' 2 ) \p x ,p 2 ) (p[ lP ' 2 \ (2.3) 
where the (4 x 4) matrix S aier2 y ia ' 2 (pi,P2,p'i,P f 2 ) is given by 

Sa 1 a 2 ,a' 1 cr' 2 (Pl,P2,p' 1 ,P 2 ) = <An<7 3 {Pi , Vi)^'^ (Pi » Pi) ( 2 - 4 ) 

Since we are only interested in the entanglement of the spin degrees of freedom we have to trace over the momenta 
of the two fermions, that is to say 

p frcc = / d 3 q 1 d 3 q 2 (q 1 ,q 2 \*)(*\q 1 ,q 2 ) 



x We could regard arbitrary spinorial or tensorial structures by simply generalizing the function ip(pi,p2, <ti, 0-2). There is no loss of 
generality with the choice made here, as we will see later. 

2 NB: In order to simplify notation we write the state \pi,P2) instead of \pi,P2) 
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as (<?i,32 bi,P2> = (0| b(qi)b(q<2)tf {pijtf (p 2 ) |0> = S(qi - ^2)^(92 - Pij^\ we obtain 

fftSiriai = J d 3 Pl d 3 p 2 S l7ia2 ^ il7/2 (p 1 ,p2,Pi,P2) (2.5) 
If we want to calculate the degree of entanglement of this kind of states we could apply any of the common entan- 



glement measures on the density matrix (2.51 



3 Asymptotic Dynamics 

3.1 Gauge Invariance and Asymptotic Evolution 

One formulation of the QED lagrangian which is very useful in order to discuss its gauge invariance is 



C = \f^F^ + i&f (fy - ieA„) i> - + ^aB 2 + Bd^A* 



Where B is a gauge fixing field [22] ■ The fields that appear in the lagrangian are (in the interaction picture) 

1>(x) = -V / d\[^J2 [bn(p)u n (p)e-^ x + dl(p)v n (p)e^} 
(2tt)2 J Mpo^ 1 



(20 i 
1 



(2n)i 

1 f rl 3 k r 
(2tt)2 J (2k )2 L 



d 3 k 



B(x) = (2n)-* / [ c \ky k -* + c{k) 



c. 



ik-x I 



/2k 

The gauge fixing field allows us to identify physical states and observables: 

• c(k) \4>) — <^4> |0) could be a physical state (since it is gauge invariant) 

• Let O be an operator corresponding to an observable, then [B, O] = 

It can be shown [2U] that in the Feynman Gauge (a = 1 ; B = —d^A^) and for asymptotic times, the field modes of 
the B field that annihilates physical states are given by the modes of the free vector potential: 

c{k) = ik ■ a(k) 

that immediately leads us to the conclusion that the states b a (jp) |0) are physical in the sense of gauge invariance. 

However, given that the large t limit of the QED Hamiltonian is not the free Hamiltonian, the states |0) do 
not correspond to the asymptotic states of QED. A well-defined physical asymptotic state must not only be gauge 
invariant but evolve with the proper Hamiltonian corresponding to the theory which describes its dynamics. We can 
conclude that the standard interaction picture does not serve to describe the correct asymptotic operators. 

Thereby, the state b\{jt) |0), with b\(p) in the interaction picture, is not adequate since (although it is gauge 
invariant) it does not correspond to a state that evolves properly at asymptotic times. 

Our goal is to find the charged particle creation and annihilation operators with well-defined momentum and spin 
that behave adequately at large t (evolving undergoing an asymptotic interaction corresponding to asymptotic QED) 
and, furthermore, that they be gauge invariant (they must commute with B-field). To accomplish this objective we 
will need to compute the interaction operator in the large t regime. 



3 We are supposing distinguishable fermions to shorten notation 
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3.2 The Asymptotic Interaction Operator in QED 



Let us show what we have announced before, that is to say, the asymptotic Hamiltonian of QED is not H , to do 
this, we begin from the interaction operator between the charged fermions and the electromagnetic field 



H int = / jnfflA fflfa = - e / : ^(x^^x) : A^(x)d 3 x 



(3.1) 



After substituting the fields expresions in (3.1l and integrate over x, the resulting expression for H lnt (t) is an 
integral over the momenta p, q, k of the fermions and the photons, whose are related by the conservation equation 
q = p+k. 



We divide the resulting terms after normal ordering expression (3.1 1 in two different groups: 



• 1st group: terms that include two creators or two annihilators of charged particle S or bd, these terms are 
proportional to the following exponential factor 



exp(±ikot) exp(ipot) exp(iqot) = exp 



i ( \J p 2 + m 2 + \j (p+ k) 2 + m 2 ± fco^ t 



2nd group: terms that include one creator and one annihilator of charged particle Wb or cftd, these terms are 
proportional to the following exponential factor 



exp(±ikot) exp(ipot) exp(-iqot) = exp 



i (^Jp 2 + m 2 — \j (p+ k) 2 + m 2 ± k ^j t 



The former acquire a phase that is highly oscillating at the limit t — > oo, thus they cancel (in the weak sense) and 
can be neglected. The latter, on the other hand, have a non-zero contribution at the limit led by the value of the 
stationary phase, which is to say, at the proximity of k = 0. 

The leading contributions to the Interaction Hamiltonian in the asymptotic regime are 



H 1 



(2w)i 



[ u m{p)l^u n (p)tf m (p)b n (p) - V m {q)^v n {q)dl{q)d m {q)\ 



al(~k) + o„(jf)] 5(p- q + k)e i ^°- p °-^ t 



d 3 pd 3 qd 3 k 
v / 2q p Q uj 



We will carry out the following approximation: let us consider k — > in all the functions that grow slowly enough 
with t, that is to say, the Dirac's Delta 5(p — q+k) will fixp=g for the spinors and the fermionic annihilation and 
creation operators, Thus we have that 



p> 



iUnip+k^Umip) —> iu n (p)l^Um(p) = $nm iv n {P '+ k)Y 'V rn (p) —> W n (p)^V m (p) = S nm 

m m 



so that we have 

nnt -e fd 3 pd 3 kp' J - 



(27T)! 



'2u> Po 



£ [ b n(p)bW) + d n (p)dl(p)} k(-fe) + o„(£) 



(p-\-k) 2 +m 2 — \J p 2 +rn 2 —u>^jt 



We could expand in power series of k the exponential and keep the leading order (as the integral is dominated by 
k^O) 



\J (p + k) 2 + m 2 — \j p 2 + m 2 — u> = \J (p + k) 2 + m 2 



k ■ V £ \/{p+ k) 2 + m 2 
k-p 



-^p 2 + m 2 - u + 0(uj 2 ) = ^ - lj + 0(uj 2 ) 

Po 



So that we can write the asymptotic interaction Hamiltonian the following way: 



d 3 k 



(2k ) 2 J (2tt)S Po 



1 d 3 p „ i£Et-iut r + , A /A 

Ve p o p{p) a T(-fc)+a M (fc) 



(3.2) 
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Which can be also written as 



if 

S+)t„\ , „(-) 



d s k 



1 d 3 p 

-v^e p o 

(27T) 3 PO P 



Where we have used that a^x) = (x) + (x) with [23] 

1 



(2tt) 3 / 2 



d*x^/2k e~ lk - x a^\x) aUk) 



(2tt) 3 / 2 



The integral over d 3 k yields to a Delta, resulting that 



where 



being 



HZ t (t) = -e J J^(x,t)a^x) d 3 x 



Po 



Po 



d 3 xy^e ik - x a,(+\x) 



P(p) = ^2 [ b n(p) b n(p) ~ dl(p)dn(p)] = p-(f) ~ P+(p) 



(3.3) 



(3.4) 
(3.5) 



In order to clarify future calculations, is convenient to write the interaction Hamiltonian (3.2 1 in the next form 

d 3 k 1 



(2fc )5 (2tt)? 



- 7 M 



aU-k) + a M (fc) 



Where 



JL= I d 3 p P -e <>(p) 
Pa 



(3.6) 



(3.7) 



Summarizing : The Interaction operator is non-zero in the asymptotic regime. The residual interaction is equiv- 
alent to the interaction between a spinless charged particle current and the electromagnetic field. The resulting 
expression is universal, and can be readily generalized considering charges with an arbitrary spin. 

Eventually, the operator J£ s has a very clear physical meaning: a given state of (free) charged particles with a 
well-defined momentum 

|*(pisi, . . . ,p n s n ,q 1 i 1 , . . .,q n i n )) = 6^ (pi ) . ..^(p^dj^gi ) . ..d\ n (q n ) |0) (3.8) 
is an eigenstate of this operator, being the corresponding eigenvalue 

m n k 

A = X^( fc ' <; *) + X^( fc ' i;p .?') where j^(k, t; p) = e— e _i ^* 
3=1 3=1 Po 

Furthermore, it is important tho emphasize that current operators in different points of the space-time commute: 

[J&(x),JH(v)] = o- 



3.3 The Asymptotic Interaction Picture 

Perturbation Theory is based on the fact of consider that, at asymptotic times, the interaction is deactivated and, 
therefore, H(t — > oo) = H thus, the "in" y "out" states that would evolve by means of the free Hamiltonian. This 
fact allow us to include this asymptotic temporal dependence within the operators, building the so-called standard 
interaction picture in terms of the Heisenberg picture: 

Oj{t, t ) = T exp (-i jf dr [H(r) - H (r)}^ H {t) T exp (i jf dr [H{t) - H (t)}^ . 

This transformation requires the introduction of an arbitrary time t in which the fields in both pictures are equal: 
Off (to) = 0/(^0; *o)- This time t is completely arbitrary since if we transform back our operator into the Heisenberg 
picture, every explicit reference to to will vanish. 
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However, we have seen that, in QED, the asymptotic limit of the Hamiltonian is not the free Hamiltonian but 
H as = fljjjt + Hq, Thus, it is necessary to build a new interaction picture in which the operators include the proper 
temporal dependence at asymptotic times. We build it as follows 

Of (t, t ) = T exp (-i jf dr [H(r) - ^"(t)]) H (t) T exp (i jf dr [If (r) - ff as (r)]^ 

Or, otherwise 

Of(Mo) = Texp (-t jf [H int (r) - H&(t)]) H (t)Texp (i ^ dr [H int (r) - ffgS(r)]) 
In order to study the dependence on to of the previous expression, let us expand the exponentials to first order 

Of(Mo) = H {t) +i [ dr (r) - # int (r), H (t)] 

J t 

Here, we clearly see that there is an explicit dependence on to that comes from the lower integration limit. This 
dependence is governed by the difference 

#£t(*o) ~ H int (t ) 

So if we choose one to that be itself an asymptotic time (large enough) the dependence on to will vanish by construction. 
Thus, we can set to = ±00. The analogous can be demonstrated for every order in the expansion. 



3.4 Building the Asymptotic Fields 

In this section, we will go from the standard interaction picture to the new asymptotic interaction picture (via the 
usual time ordered exponential of the difference between the complete Hamiltonian and the asymptotic one) Thus, 
the field operators will include the asymptotic evolution. Given the functional form of the asymptotic Hamiltonian, 
the fact that the asymptotic current operators commute and that [a' , a] is a C-number, it can be proved that 



T exp (i jf drH^ir)^ = exp jf drH^r^j exp Q jf d n £ dr 2 [H&(n) , H&fo)] 
We can now obtain the field operators (fermionic and photonic) in the asymptotic interaction picture. 



(3.9) 



Concerning the photon field, provided that [[H^ t (ri), H^fa)}, AJx)\ — 0, is easy to build the asymptotic vector 
potential operator, 



so we have that 



AZ(x)=exp(iJ t drH^(r)^ A^x) exp (-ijf drH&( 
A» s (x) = A M (x) - e f drd 3 yD(r - t,y- x)J^(r,y) 

Jtn 



It is important to remark that, from this expression, we obtain DA^(x) = —eJjf(x), that is to say, this asymptotic 
vector potential contains the Coulomb field generated by the asymptotic charged current Jjf plus the free field. This 
is reasonable since at asymptotic times the residual interaction is like the interaction between a charged current and 
the electromagnetic field. 

Furthermore, given that J^ix) commute with itself, the asymptotic bosonic field obeys the same commutator as 
the free field 

[A™(x),A™(y)] = -i giiV D(x-y) 



Let us consider the matter field, in this case both terms in (3.9 1 contribute. The computations are not difficult but 
a little bit tiresome (above all for the phase factor) , the detail can be found in [T^] and . The result is that there 
exists a transformation U(t) that goes from the standard interaction picture to the asymptotic interaction picture 
which is given by 

U{t) = exp [R(t)j exp [i$(t)] 
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Where we will call the operators exp [R{t)] and exp [i$(t)] distortion and phase operators respective!}^] which are 
given by 



R(t) 



(27T) 



$(*) 



aUk)e "o f 



pk 



a M (fe)e 

P-Q 
yj (pq) 2 - m 4 



p(p)dp 



dk 



(2k )i 



sign(t) In — dpdq 
to 



(3.10) 
(3.11) 



3.5 Asymptotic Evolution of the two Particles State 



For our purposes, we would like to define a two charged particles state (eigenstate of the fermionic number operator 
with eigenvalue 2) with well-defined momenta and spin projection. In the standard interaction picture such a state 
would have the form 

|*) = 6i 4 (pa)&kCPi)|0) 
where the operators include the free evolution within. 

Nonetheless, as we argued above, the interaction picture is not adequate to describe the dynamical behaviour of 
charged particles in the large t regime. If we transform the operators to the asymptotic interaction picture we will 
have 

|*) = bf 2 Hp 2 )b^( P i) |0) = exp [R(t)] exp [»*(*)] ^(paj^fa) |0) 

Taking into account that the fermionic Fock states are eigenstates of the asymptotic current operator, It can be 
obtained an expression in which the asymptotic factors are free from fermionic operators 



j<Ku r ,t) e w( Pl , P2 ,t) b ^ (pi)6 t (p2) | Q) 



Where 

W( Pl , P2 ,t) 



(27T) 



d 3 k 



kp 1 



pik 



kp 2 



U 

Pik 



aUk) 



pik 



p^ 
Pik 1 



a^k) 



47T 



u r 1 (pi,p 2 )log 



to 



(3.12) 

(3.13) 
(3.14) 



are the "eigenvalues" of the R(t) y Q(t) operators associated to the considered state [TB] and where u r (pi,p 2 ) is the 
absolute value of the relative velocity between the two particles, u(p, q) = yl — m 4: /(pq) 2 . 

It is obvious that this construction is not Gauge invariant since W(pi,p2,t) contains exponentials of photon 
operators. This fact could make us think that the charged particle concept could not be well-defined. 



Furthermore, the resulting phase factor (3.141 (that is not often considered in the literature as it is of second order 
and it does not appear in the cross sections) Although it is gauge invariant, it seems to be divergent for to — ► ±oo. 
In the following sections we will demonstrate that this apparent divergence disappears when we apply the "dressing" 
procedure building physical gauge invariant states. 

To end this section let us mention an interesting result that we will use later: It can be proved |16j that the 
transformation of the asymptotic state iV'as) ~ > IV'D such that 



|V4) = e *K.') e ^w.<*,«*,t) 6 t ( Pl )bUp 2 ) |o) 



(3.15) 



where 



W(p 1 ,p 2 ,c 1 ,c 2 ,t) 



(2n)i 



11 

Pi 
Pik 



d 3 k 

J 1 



Pi 
Pik 



„1 


+ {pik 






&-*} 


e P2 o 







aUk)- 



a^{k) 



4 Notice that we have set to — > ±oo through the calculations due to the arguments argued in section [3^3] (as it is directly done in |16| . 
|20| and [23]) except for the phase operator, in which we will keep the explicit dependence on to because its contribution is apparently 
divergent for asymptotic values of to- We will show that the dependence on to in the phase operator is irrelevant, and the dressing 
procedure will cancel it. 
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(that differentiates from W in the terms and Cj) operates within the asymptotic states Hilbert 7i a s, that is to say 

l^as) G U &s & |Va S > G ^as 

If the following integrals converge 



— e p o 
pk 



2ko 



and 



qk 



e "o 



(3.16) 

This is a formal result which we will use in this work to understand certain aspects of the "dressing" procedure. 



4 Building Physical Asymptotic States 
4.1 The "dressing" procedure 

The states built so far (that evolve correctly in the asymptotic regime) are not physical (they are not gauge invariant). 
The solution to recover the physical particle conception [20], [21] consists on "dressing" the fields by means of operators 
that restore the gauge invariance without altering their dynamics. This procedure that allows us to define a physical 
particle state in the asymptotic regime. 

The idea is to "dress" the matter field with a "dressing" operator that under gauge transformations behaves 
opposite to the field and hence the built object will be gauge invariant. This is to say that given a field <p(x) that 
under gauge transformations behaves as follows 

<p{x) — ► ip{x)e le9{x) 

we need a dressing such that 

h-\x) —> h- 1 (x)e- ie0 ^ (4.1) 

so that the dressed field 

= h-\x)(p(x) (4.2) 

is gauge invariant. 

But this is not the only condition that the dressing should satisfy: We will also demand that it conserves the 
dynamics of the field that we are dressing. First of all let us consider a free matter field in the infinite-mass limit it 
can be proved [24] that velocity is superselected and the equation of motion for the field has the universal form 

u ■ d<f>(x) = (4.3) 

where u is the four- velocity of the heavy particle. This equation is simply the statement that the field is constant 
along the world line of a particle movin with 4- velocity u. Consequently, if we parametrise the world line of a particle 
that moves with a 4-velocity = 7(77 + v) where 

• r\ is an unitary temporal vector 77 = (1,0) 

• v = (0, v) being v the 3-velocity of the particle 

2\-V2 



• 7 = (1 — \V 

we obtain that 



x"(s) =x"+ ( a - x°)(ti + v)" (4.4) 



Then, (4.3) implies that, for an arbitrary s, $ [x(s)] = $ [#(0)]. 

If the matter field is minimally coupled to the EM field the equation of motion turns into 

u ■ D$(x) = where = <9 M - ieA^. (4.5) 

This analysis, which is done for the heavy matter sector, is also valid for any field that asymptotically could 
correspond to a charged particle with four-velocity u. We will give two arguments to justify this statement: 
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1. The asymptotic dynamics of QED is governed by soft photons for whom any electron is heavy. 

2. It can be shown [25J that the asymptotic interaction Hamiltonian vanishes for the propagator of the dressed 



fields that satisfy (4.5 1 if one is at the correct point in the mass-shell. 



If we demand the dressing preserves the heavy particle dynamics (we impose (4.51 to (4.2 1) we will obtain the next 
equation for the dressing operator [5D] (that is called "the dressing equation"): 

u -dh^ix) = -ieh- 1 (x)u- A(x) (4.6) 



Summarizing, the equations (4.1) and (4.6 1 are the two fundamental requirements that the dressed field $ must 



satisfy in order to describe physical charged particles. 
4.2 "Dressing" the two particles state 

The dressing operator hT 1 consists of a phase term and a distortion term (so the asymptotic evolution operator does) 
[2D] that is to say, we can express 

h- x (x) = e Xi ^e- iKi{ - x) 
where it can be shown [20] that in the large t limit 

r , [ d 3 k 1 / _ ik . x Walik) lk .\ 

where 

V? = to + v i y{r)-v i )'k-k> 1 

And also that 

e- iA '' =exp{-ie£(r 1 + v i y d ^[x( S )]ds} (4.8) 

where we have an integral along the world line of the massive particle with constant 4-velocity parametrised by 
means of s [20] 

The action of the operator -^q is defined [20] , [21] as 

1 -J(x) = I d 3 zG(x-z)f(z) where G(x) = - [ ,f* } - ' 



Q-d } J (27r) 3 fe2- (ff.fc)2 4tt y/S^ + ^v-x) 2 ' 

Given that the dressed matter field has the expression <$>(x) = h~ 1 (x)(p(x), the dressed creation operators in 
function of the dressed asymptotic field is given by 

&Ld(P> ff )= / d 3 x.[^ui(p)e- lli '^e^ x ^ aa (x)e^ x 
J V Po v v ' 

that is to say 

Expressing the asymptotic creation operator in terms of the creation operator in the standard interaction picture: 

bljq,*) = (2-)- 3 E / A' J^uMuA't) J ^e^''^'W e ^'W ge^e^ 

bl(q') 

It can be demonstrated that e - lK i( x ) commutes with e^^e 1 ^^ (since current operators commute) so that 



&L,dM = ( 27r )~ 3 £ / d 3 q'J^u a (q)u a ,( q ') J d 3 xe^^e^e R We^e-^bHq\a')e"^e- R ^ 
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Let us start computing the phase operator: The relevant dynamics of the gauge fields in the asymptotic regime 
is given by [25], [22]: 

where we see that the matter coupling comes from the asymptotic current. Using this and omitting the dependences 
with the gauge fixing field, that are irrelevant acting on physical states, estados ffsicos, we can write Ki as 



^2 ft 



Ki 



4tv 



ds(f] + ViY \ d 3 z 



y/{x(s) - z} 2 + 7 2 [w ■ {x(s) - z}]< 



Substituting the asymptotic current (3.5 1 we obtain that 

„2 ft 



Ki = 



-Itt 



to 



ds(r] + ViY / d 3 pj—p(p)f(s,v i7 Xi) 
Po 



where 



f(s,Vi,Xi,Pj) 



x(s) — s 



pj_ 



Vi ■ x(s) + s 



Pt 

Pjo 



21 ~ 



(4.9) 



(4.10) 



We need to commute e iKi ( x > with b\ in order to do it, we will use Hadamard's Lemma: 

e xA Be - xA = B+ [A,B]x+j [ [A, [A, B]]x 2 + ... 



Since we have that 



Let us use Hadamard's Lemma with 



-iK, 



b\q>) = e- iK ^(q') 



e 'e 



A = Ki B = tf{q') 



>^tf(q')e iK * = (^{ q ')-i[ Ki ,b\q')\+{-if^ [K^tf)]] + 



And the commutator is 



Since we know that 



.2 ft 



[Ki,bHq')] = / {fi + ViY / d 3 P1 ^ [ P (p),bHq')] 
J to J Po 



[p(p),b\q')] = [tf(p)b(p),tf(q')]=tf(q')6(p-q') 
we can perform the integration over d 3 p, resulting 

= 2 ft 



[Ki,b\q')] 



-in 



(r) + ViYj^tf(q')f(s, v^x^ds = k(vu q', 



(4.11) 



On the mass shell of q', it is also satisfied that Vi — (f /q' , so that, (given that q^/qo — (»/+ v) M and (rj+v^^^+Vi)^ 



e 2 1 



dsf(s,v i ,S i ,q) 



4tt 7 J to 

Therefore, the dressed asymptotic creation operator has the form 



^ > d(9^) = (27r)- 3 E/ d V 1 



E' q E q 



(q)u a ,(q') / d 3 a;e < («-«'>e x «W e fl(*) e »(t) e -f«*W 6 t(g/ j£r /) e -«f«(») e -fl(*) e -»(t) 



For our purposes we need to calculate the product of two dressed asymptotic creation operators applied to the vacuum 
state 



&L,d(9l^l)&L,d(9l^l) 10) = (27rr 6 ]T / d 3 q[d% J m ' ui(q)uAq'H(q)uAq') ' 

CTCr ' J \ 1l Pi 12 P2 

d 3 Xl d 3 x 2 e l{qi ~ q 'i )xi e %{q2 ~ q i )x2 e Xl{x) e X2{x) e R{t) e^e |0) 
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To arrive to this expression we have used that e~ iKi |0) = |0), i = 1,2 and we have performed two commutations^] 
Now, we need to commute e~ lKl ^ with b^(q 2 ), the computation is analogous to the other two performed commu- 
tations with the only difference that the velocity that appears in the operator is on the mass shell of q[ and we are 
commuting with the creator with momentum q' 2l the calculations yields to the term 



Kl2 = |-7^ ds(ri + vi)^(r) + v 2 )^f(s,vi,x 1 ,q' 2 ) = |~7(] r\-r 2 ) j ds ./'(■-. n. . ,r, . </ 2 ) 



Therefore, the two charged fermions state has the form 



&L,d(9i^i)C,d('?i^i)|0>-(27r)- 6 ^ f d 3 q[d% J m ' * uUq)u a ,(q>Uq)uAq')- 

aa i •* \ 11 Pi 92 P2 

d 3 Xl d 3 x 2 e l{qi - q 'i )xi e l{q2 - q 'i )x2 e Xl{x) e X2 ^ |0) 



The next step is to compute the integrals over the affine parameter s along the world line of the massive particles, 
in order to compute the phase factors we will parametrise it in the following way 



x»(s) =x fi + (s-t)(r) + vY 



x(s) = x + (s — t)v 
x°(s) = s 



4tt 7 J tQ I J 



integrating trivially 



e 2 1 



t-t 



An 7 yj(x - v'tf + 7 2 [v> ■ (x - v't)f 



The term in the denominator is a retarded position R = x — vt analogous to the retarded position in the Lienard- 
Wiechert potentials HHI (we will discuss this point later): 



e z 1 



t-t 



(4.12) 



At large times, the leading contribution to R is the term —vt since we are considering the region of the space- 
time in which we have large times but finite positions as in the spatial infinite the interaction actually vanishes (see 
discussion in [20 ). Taking it into account, the contribution of x is subdominant and we can assume that we are 
always at positions given by the asymptotic times so that the phase can be factored out from the integral over x in 
the asymptotic regime approximation. 



Analogously, we obtain that k,i is 



«2 



e 2 1 

47T 7 



t-tn 



Rl + 1 2 & • R 



(4.13) 



Let us show that the phase factors (4.12 ), (4.13 1 have a clear physical interpretation, writing the Licnard-Wiechert 



potential created by a moving charge with 4- velocity u^ L = (7, u) at large distances (cf. [29] Cap. 14): 



.4 



class 



4tt 



R 2 



u-R 



-in 



i? 2 + 7 2 



v-R 



Where R is the retarded position. On the other hand, the asymptotic current is 

J* s (t, x) = J d 3 P (n + v)^p{ P )5 3 (x - vt) 

5 We have commuted e~ xKl with b' (q^) and e~' K2 with 6' (q'2), giving us the two C-numbers e~ lK1 and e~ 
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acting with it on a state b'(p) |0) we obtain the next eigenvalue: 

j™(t,x) = (n + v)^ 3 (x-vt) 
Thus, the effective coupling classical field-asymptotic current acting on a state b^(p) |0) results 

S 3 (i?) 



e 1 
47T 7 



i? 2 + 7 2 



v-R 



Where has been used that (n + v)^(n + v)^ = I/7 2 , Therefore, we can write 



That is to say, These phase factors can be identified as the temporal evolution due to a Hamiltonian that comes from 
the coupling of the asymptotic current with one classical Lienard-Wiechert potential created by a current at a large 
retarded position (large at asymptotic times) and governed by the term Vit. 



We have now to calculate the crossed phase term k 12 which is given by (4.141 

K12 = — 7 (1 — vi-v 2 ) I ds |[x - Hit + s v r ] 2 + 7 2 [v\ ■ (x — vit + s u r )]j 



-1/2 



Or, as a function of the retarded positions, 



Ki2 = ^7(1 -vi-v 2 ) I ds 



Ri + sv r 



vi ■ [Ri + sv r 



-1/2 



We can decompose the retarded position in two components, one normal to the relative velocity and the other in its 
same direction 

R = b + vi — vi\\ + vi± 

Thus, we will obtain the following result (operating, now, with the modulus) 



Ri + sv T 



= b 2 + (R ll + s Vr y 









2 




v 2 ■ 


(Ri- 


\-sv r j 




h ( 



Vl± 



b+Ru+ sv, 



For simplicity, let us also consider that vu_ — =>■ Hi — vi\\, situation that would be satisfied, for instance, working 
in the center of mass frame or a frame in which cither particle is at rest. 



«12 



4tt 



7(1 -Vi-v 2 ) I ds{b 2 + (R\\ +sv r ) 2 + j 2 [vi (R\\ + si> r )] 2 J 



-1/2 



performing analitically the integration we eventually obtain 



Kl2 = -r-7(l - Vi-V 2 ) 

47T "fV r 



In h In 

to 



R \\ 
t v r 



b 2 h 2 



{tv r + R\\Y 



R\\ 

t v r 



b 2 h 2 



(to v r + R\\) 



For asymptotic t, to we can neglect the quotient terms in the second logarithm according to the arguments viewed 
above in this section, resulting thai[^] 



e 2 1 t 
«12 ~ -rl (1 - vi -vs) In — 

47T JV r Iq 



6 Although Ru is dominated by —v\t it can be shown that, in any case, the asymptotic behaviour of the factor K12 is correct. To see it, 
it is only needed to set in a frame in which v\ = and take into account that x/t — > in the asymptotic regime as we have argued above. 
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On the other hand, the expression of the relativistic relative velocity between two particles as a function of the 
3- velocity of each one is 

v r 



1 - Vi-V2 



So, substituting it in the previous expression we obtain 



e 2 1 , t 
K\2 = In — 

4-7T U r to 



(4.14) 



Thus, the two particles state results 



&L,d(?l.*l)&L,d(«l>*l) 1°) = / \\ E , Er uU q )u a ,( q i)uUq)u a ,( q >) ■ 

aa' V qi P' 1 92 P 2 

. e -i Kx {vi,t) e -i K2 {v' 2 ,t) I d 3 Xi d 3 X2 e i(gi-gi)^i e i(92-^)x 3e ii(i) e xi(x) e x 2 (x) e -m 12 (a ; ) e i* 6 t(^) & t(^) |Q) 



But, recalling the section 3.5 



e i *6t(^) 6 t(^)|o) = e <*6t(^) 6 t(^)|o) 



where is given by (3.14) and is identically equal to — k 12 and, therefore, both phase factors cancel themselves, being 
our state 



C(^i)O^) i°) = ( 2 ^r 6 E / J E , m uu q )u a ,( q ')ui( q )u^( q ') ■ 

aa > J \ 11 Pl 92 P'l 

. e - lK i(Rut) e - lK2 (R 2 ,t) e w(t) f d3 a . 1 dS a . 2 e*(« 1 -«i)«i e i(«i-^)*>eXi(*)e xa W6 t (^ 1 )6 t (^) 



where we have already applied e R ^ on the two particle state, and W(t) (given by (3.13)) has no fermionic operators. 
In order to arrive here we have used that e R ^ and e Xi ^ x ' commute (It can be readily probed since the first commutator 
of the BCH formula vanishes) 

We have now to consider the product of the operators e Xl ^e X2 ^ x ' and to integrate over d 3 Xi and d 3 q [. 

It is shown in [5U] that the result of the integration is (at asymptotic times)is the following 

&L,d(^)&L,d(^) 1°) = e- tel(fll ' t) e- i ^(*'*)e^W^( Plj ^ ) t)/i^(p 2> ^,06 t (8i^ 1 )6t(g i , ffa ) |0) 

where h~^ {t (pi, m, t) is called the "minimal part" of the dressing and it has the following expression 



fy 3 oft(P*'*> l; *) = ex P 



d 3 k i / vr 



(2tt) 3 / 2 ^2ko \Vi-k 



a M (fc) 



V 



Vvk 



<(k) 



where 



V? = {r) + v i y{r l -v i )-k-k> 1 



Using the Baker Campbell Hausdorff relations we can demonstrate that 
where 



W(p 1 ,p 2 ,v 1 ,v 2 ,t) = 



d 3 k 



(27r)5 J VW 



and where 



vr 



Vvk 
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Notice that even the first commutator in the BCH formula vanishes in this case. 

So we can conclude that the dressing procedure is equivalent to perform a transformation on the two particles 
state that operates within the asymptotic states Hubert space since cf satisfy the conditions (3.161. 

Moreover, it is shown in [20 that operating on physical states on the mass shell, (that is to say, p^ = m'y i (i]+Vi) ,i ), 
the operator exp(W') is only function of the gauge fixing fields, thus, between physical states (which are annihilated 
by c) on the mass shell, this operator acts as the unity. We can, hence, ignore it when we are studying our case of 
two physical particles on the mass shell. 

Eventually, we have obtained a final expression for our dressed asymptotic state of two charged particles 

K) = e-^^e-^'H^bl^) |0) (4.15) 



where n{R\,t) and n(R 2 ,t) are given by the expressions (4.12) and (4.131 and they have the interpretations discussed 
above (evolution undergoing an effective interaction asymptotic current-classical field) 

The result is completely symmetric, we could repeat the computation inverting the order of the fermionic operators 
and the result would be exactly the same. 



5 Entanglement of a dressed asymptotic sate of two spin 1/2 charged 
fermions 



Let us build an arbitrary dressed asymptotic state of two charged particles analogously to (4.151 
|*)= 53 / d3 Pi d3 P2^ 1( T 2 (Pi,P 2 ) 



Where we are using the notation (2.2 1 



The density matrix associated to this state is 



&1,&2 a\ ,er' 



where 



e i8 = e -in(R u t) e -in(R2,t) e in(R[,t) e iK(S^,t) 



which does not depend on the spin of the fermions. Thus, we can write the previous expression as a 4 x 4 matrix in 
the spin space 

|*><*|= f d 3 Vl d 3 p 2 d 3 p\d 3 p' 2 ^ S^y^p^p^p'^ (5.1) 



Being S aia2 y i<J ' 2 (p ll p 2l p' 11 p' 2 ) the matrix ([2^ I. 

Since we will only consider spin entanglement, we have to trace over the fermions momenta 



J d 3 qi d 3 q 2 ( qi ,q 2 |tf ) (y\ qi ,q 2 ) 



Using that essentially (qi , q 2 \pi , p 2 ) = (0| b{q{)b{q 2 )tf (pi)tf (p 2 ) |0) = 5(qi — p 2 )S(q 2 —pi), those deltas fix p = p' and, 
hence, given that Ri ~ — V{t, they fix — Ri. Therefore, in the asymptotic limit, e lB = 1 and all the phase terms 
vanish: 

Pspin = J d 3 p 1 d 3 p 2 S ai<T2i<J > i(T i 2 {pi , p 2 , pi , P2) 



this expression is equal to (2.5) 



Thus, this is the same result which we would have obtained without considering any effect of the infrared structure, 
so, eventually, we have obtained that the spin entanglement is not modified by the effects derived from the presence 
of soft-photons 
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6 Conclusions 

We have calculated the effect of the infrarred structure of the QED on the spin entanglement. 

In p] , it is pointed out the ignorance of the effect that the unavoidable presence of an undeterminate number of 
soft-photons could have on the entanglement of charged qubits, questioning if the very concept of qubit is no other 
thing that an idealization. 

In order to approach this problem we have used the dressing formalism to, first of all, build physical states of two 
charged particles from the proper asymptotic states of QED and, afterwards, analyse what happen with their spin 
entanglement. 

As the main result of this work, it has been demonstrated that the infrared structure of QED has no effect on the 
spin entanglement of charged qubits, thus all the spin correlation functions (and in particular the EPR correlations) 
are not modified by the presence of soft-photons when we are regarding as qubits the proper physical asymptotic 
states. 

Finally, Along this work the next contributions have been made 

• It has been discussed the physical relevance of the parameter to that appears in the building of the new 
interaction picture. 

• It has been shown that the dressing procedure operates within the asymptotic states Hilbert that is to say 

h~ Has C Has- 

• The phase contributions of the dressing procedure have been explicitly calculated for the first time. 

• It has been shown how the dressing procedure cancels the unphysical critical phase dependence on to that comes 
from the construction of the asymptotic states. 

• The effect of considering the infrarred structure of the QED on a two charged fermion state has been calculated, 
and it has been shown that the additional terms that appear are equivalent to the temporal evolution undergoing 
an effective interaction Hamiltonian between the asymptotic current and the classical EM-field. 

In the future we will investigate what are the effects of the infrared structure of QED regarding the momentum 
entanglement or the cross correlation spin-momentum and, summarizing, we will try to determine the effects of the 
infrarred structure of QED on all the quantum information tasks. 
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